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Abstract. A finite stack of thin superconducting tapes, all carrying a fixed current 
/, can be approximated by an anisotropic superconducting bar with critical current 
density Jc — Ic/^aD, where Ic is the critical current of each tape, 2a is the tape 
width, and D is the tape-to-tape periodicity. The current density J must obey the 
constraint / Jdx = I/D, where the tapes lie parallel to the x axis and are stacked 
along the z axis. We suppose that Jc is independent of field (Bean approximation) 
and look for a solution to the critical state for arbitrary height 2b of the stack. For 
c < \x\ < a we have J = Jc, and for < c the critical state requires that = Q- 
We show that this implies dJ jdx = in the central region. Setting c as a constant 
(independent of z) results in field profiles remarkably close to the desired one (Bz = 
for |a::| < c) as long as the aspect ratio hja is not too small. We evaluate various criteria 
for choosing c, and we show that the calculated hysteretic losses depend only weakly 
on how c is chosen. We argue that for small D j a the anisotropic homogeneous-medium 
approximation gives a reasonably accurate estimate of the ac losses in a finite Z stack. 
The results for a Z stack can be used to calculate the transport losses in a pancake 
coil wound with superconducting tape. 

PACS numbers: 74.25. Sv,74.78.Bz,74.25.0p,74.25.Nf 

Submitted to: Institute of Physics Publishing 
Supercond. Sci. Technol. 



§ To whom correspondence should be addressed. 



ac losses in a finite Z stack 2 



d 

D 



k 2a >[ 

I I 

Figure 1. Finite Z stack: a stack of superconducting strips of infinite length in the y 
direction, each carrying current /. The overall height of the stack is 2b. 

1. Introduction 

It is now possible to wind pancake coils from long lengths of high-temperature 
superconducting (HTS) tape P,|2]. To date there have been no theoretical calculations 
of the losses in this difficult geometry other than those using a variational approach [3] 
or via numerical simulations [2]. A related geometry, shown in figure 1, consists of a 
stack of tapes of infinite length in the y direction, each carrying a total current /. This 
is closely related to the coil geometry but has the advantage of computational simplicity. 
The tapes of width 2a in the x direction are stacked in the z direction to a height 2b. 
This problem should be distinguished from one that has been previously considered [1], 
where the total current carried by all the tapes was given but the share taken by each 
one could vary, as would be the case if the tapes were bundled together to increase 
the current- carrying capacity of a conductor. In a coil the current in each winding is 
constrained to be the same, and we preserve this feature in our related geometry. We 
assume further that the thickness d of each superconducting film is much smaller than 
the width 2a. This is certainly true for 2G (second generation) YBCO tape [5], where 
the tape-to-tape periodicity D is governed by the thicknesses of the insulator, substrate, 
and buffer layers. 

Both Mawatari [6] and Miiller [71 [sjj jl considered this problem in the limit 6 — oo 
and obtained analytic expressions for the fields and currents. Both authors noted 
that in the limit D <^ a the solutions approach those for a uniform infinite slab 
of width 2a carrying an average current density I/2aD. In other words, the stack 
becomes equivalent to a homogenous superconducting slab with critical current density 

11 Reference [5] also corrects typographical errors in expressions for the ac losses given in [7]. 
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Figure 2. ac losses in an infinite Z stack, calculated from the analytic solutions in [6], 
normalized to the losses in the equivalent uniform slab, at various current amplitudes 
/ and stack periodicities D. 

Jc = -^c/2a-D, where Jc is the critical current in each tape. We expect that in practical 
applications the ratio D/a will lie in the range 0.01-0.2. In figure 2 we show exact 
calculations of the ac losses as in [6], normalized to the ac losses calculated using the 
homogeneous approximation, as a function of D/a. It can be seen that the homogeneous 
approximation is reasonably accurate for small D/a and large I/Ic- Specifically, if 
we restrict ourselves to currents of amplitude greater than 0.2 J^, this approximation 
gives better than 20% accuracy if D/a < 0.2. From an engineering perspective this 
sort of accuracy is usually adequate, especially since the error is in the right direction 
(overestimating, rather than underestimating, the dissipation). 

At present there are no analytic solutions available for the problem of a finite stack 
of conductors. To initially approach this problem it makes sense to use an approach 
that has some of the features of a homogeneous model. However, our model must 
also account, at least approximately, for the screening by subcritical portions of the 
superconducting strips. It is likely that the error in this approach will be similar to 
that of the infinite stack; see figure 2. The current density Jy and magnetic induction 
B are averaged over a volume D^; that is, we use only macroscopic values of these 
quantities. To model the constraint of constant total current in each tape, we require 
that / Jydx = I/D for all \z\ < b. In section 2 we use this anisotropic homogeneous- 
medium approximation to calculate the ac losses of a finite Z stack of superconducting 
tapes. We discuss and summarize our results in section 3. 
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2. Anisotropic homogeneous-medium approximation 

We consider a sample initially in the virgin (magnetic-flux-free) state and examine the 
initial penetration of magnetic flux as current is applied in the y direction. We anticipate 
that, similar to case of an inflnite slab, we will have a region c < \x\ < a with Jy = J^- 
For simplicity, we use the Bean P, [10] critical state model, in which is independent 
of fleld. Unlike the behavior in a homogeneous inflnite slab, however, in principle we 
should allow for c to vary as a function of z. Further, we cannot assume that = 
and B = in the region |a;| < c(z), as is the case for the homogeneous inflnite slab. It 
is known from studies of the critical state model in an isolated superconducting strip 
p!T|[T2] that no signiflcant amount of magnetic flux can penetrate subcritical portions of 
the strip (i.e., = wherever Jy < Jc); this is also true for each of the strips in the Z 
stack. On the other hand, a flnite is allowed, since magnetic flux can thread between 
the superconducting layers from the ends of the tapes without fully penetrating any 
superconductor. This leads to important constraints on B and Jy = Jm in the middle 
region |x| < c{z): Since V ■ B = 0, we must have dB^/dx = 0, such that B^ depends 
only on z. Ampere's law requires that fioJm = dB^/dz — dBz/dx. Since the second 
term on the right-hand side is zero and the flrst depends only on z, we conclude that Jm 
can depend only on z. Thus the current density Jy as a function of x has a step-function 
character, with the values Jm for |x| < c{z) and Jc for > c{z). To have a flxed total 
current in each layer we require 

JmlJc = l-{a/c){l-Ilh). (1) 

For flnite values of 6, the current density Jc in the region c < x < a contributes, via 
the Biot-Savart law, a positive value of i?2(c, 0), while the current density Jc in the 
region —a < x < — c contributes a negative value of smaller magnitude. In order to 
make ^^(c, 0) = 0, the current density in the region — c < x < c must obey Jm > 0, 
so that it makes a negative contribution to Bz{c, 0), thereby cancelling the net positive 
contribution from the currents in the regions for which c < |a;| < a. Since < Jm < Jc 
and < c < a, we thus see that c/a can vary in the range from (1 — ///c) to 1. In the 
limit as 6 ^ oo, we must flnd that Jy/ Jc — > for |x| < c and that c/a ^ 1 — I/Ic- 
The theoretical problem thus reduces to flnding a c{z) that yields macroscopic flelds 
consistent with the above requirements of the critical state. This means that we must 
have a region deflned by |x| < c{z) where Bz = 0. 

Our primary goal in this paper is to calculate the hysteretic ac losses in a Z stack. 
Using the above approach, once we obtain the solutions for Bz{x, z), we begin by flnding 
Q'inity the energy per unit length dissipated upon initial penetration of magnetic flux, 
i.e., when the current in each tape is raised from zero to a maximum value / < Jc, 
starting from the virgin state (no trapped magnetic flelds in the superconductor). To 
derive Q'^niv (^) integrate J • E over the cross section of the stack, (b) neglect the 
relatively small losses in the tapes in the middle region, |x| < c{z), where Jy < Jc, 
Bz{x, z) = 0, and E = 0, (c) note that Jy = Jc in the outer regions, c{z) < \x\ < a, (d) 
apply Faraday's law, V x E = — 9B / dt, (e) integrate over time as the magnetic induction 
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B increases from zero and reaches its final value, (f ) make use of the symmetry that the 
losses are the same in all four quadrants of the xz plane, and (g) do a partial integration 
over X. The result is [31 [TJ [131 E] 

pb ra 

Qinit = ~^Jc dz dx{a — x)Bz{x, z). (2) 

JO Jc{z) 

The physical interpretation of this formula is that Q'inn is the summation of the energy 
dissipated by vortices as they move a distance a—x from the edge to their final positions; 
the force per unit length is (poJc, where (po = h/2e is the superconducting flux quantum, 
and the density of vortices is B^^x, z) / (pQ. The current and field distributions during 
the initial penetration of magnetic flux are not the same as those that occur during one 
quarter of the ac cycle. Nevertheless, it can be shown [151 [16] that Q' , the hysteretic ac 
loss per cycle per unit length, is given by Q' = '^Q'init- 

For an infinite slab of thickness 2a, the magnetic induction upon initial penetration 
for X > c = a(l — ///c) is B^^x) = —^qJc{x — c), and the hysteretic ac loss per unit 
length associated with a cross-sectional Aah is 

= ll^oJla%I/hf. (3) 

A good starting point for the calculation of the magnetic induction inside the Z 
stack is the assumption that c = const, independent of z, and this is the approximation 
that we shall use for the remainder of this paper. Expressions for the magnetic induction 
B(a;, z) = xB^{x, z) + zB.^{x, z) = V x A and the corresponding vector potential A = 
yAy{x, z) generated by a current density in the stack Jy = = — (a/c)(l — I/Ic)]Jc 
for < c and Jy = Jc for c < \x\ < a can be obtained by using the Biot-Savart law and 
integrating over the cross section of the stack, |x| < a and \z\ < b. Results obtained for 
a constant value of c are given in Appendix A and Appendix B, and we have evaluated 
them numerically using Mathematica [17]. It can be shown that these expressions can 
never exactly satisfy the requirement that Bz{x,z) = for all |x| > c and \z\ < b. 
However, if we choose c/a to make the average of Bz{c, z) over the region < x < c and 
< z < 6 equal to zero [see flB.5|) ]. we find that Bz{x, z) ~ for all < c and \z\ < b 



with an accuracy that improves as 6 — oo. Figure 3 shows the dependence of c upon 
the current I for various stack heights 26, and figure 4 shows corresponding plots of Jm 
([I|) vs current. As expected, in the limit as 6/a — > oo, c approaches the limiting value 
a(l — I/Ic) for an infinite slab of thickness 2a and Jm approaches zero. 

Figures [5l [6l [3, and [8] show contours of constant Ay{x, a) for values of c obtained 
from ( IB. 51) . Since these contours correspond to magnetic field lines, an exact solution 
would have all these contours parallel to the x axis in the middle region |x| < c, where 
we should have By{x,z) = 0. The degree to which these contours meet this criterion 
is one measure of the accurary of our method of approximation. Despite the simplicity 
of our approximation, the magnetic field lines are remarkably straight in the middle 
region \x\ < c between the vertical dashed lines, especially for stack aspect ratios 
b/a > 1. However, as can be seen in Figs. 7 and 8, the field lines deviate from the 
desired straightness in the middle region for smaller aspect ratios. 
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Figure 3. The constant c (in units of a) obtained from (IB.SP as a function of I/Ic for 
stack aspect ratios b/a — 0.2,0.5, 1,2,5, and oo (top to bottom). For the latter case, 
c/a = 1 - I /la. 
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Figure 4. The reduced current density j„ — J,nl Jc in the middle region < c and 
|z| < 6 as a function oi I/Ic, calculated from (jB.ip and (|B.5|1 . where c is assumed to be 
independent of z, for stack aspect ratios b/a = 0.2,0.5,1,2,5, and oo. For the latter 
case, jm = 0. 
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Figure 5. Contour plot of Ay{x,z) vs x and z (in units of a) inside the Z stack for 
I/Ic = 0.5 and b/a = 2. The contours correspond to magnetic field fines flowing in 
the clockwise direction. The vertical dashed lines mark the boundaries of the middle 
region at \x\/a = c/a = 0.5895. 

In figure 9 we plot Bz{x, z) at various heights z above the center hne for 6 = a and 
I/Ic = 0.5. Note that although \Bz{x, z) \ is generally much smaller in the middle region 
< c than in the regions carrying a critical current (c < < a), it is not precisely 
equal to zero, as would be the case for an exact solution. Note also that although 
\Bz{x, z)\ is very nearly zero in the middle region < c for ^ = 0, it deviates from this 
behavior as we move away from the center line. 

Several different criteria could have been used to determine the constant c. Our 
choice, based on the best appearance of the contours of constant Ay{x,z), is to use 
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Figure 6. Contour plot of Ay{x, z) vs x and z (in units of a) inside the Z stack for 

///(, = 0.5 and b/a = 1. The contours correspond to magnetic field lines fiowing in 
the clockwise direction. The vertical dashed lines mark the boundaries of the middle 
region at \x\/a = c/a = 0.6338. 
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Figure 7. Contour plot of Ay{x, z) vs x and z (in units of a) inside the Z stack for 
I/Ic = 0.5 and b/a = 0.5. The contours correspond to magnetic field lines fiowing in 
the clockwise direction. The vertical dashed lines mark the boundaries of the middle 
region at \x\/a = c/a = 0.6809. 
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Figure 8. Contour plot of Ay{x, z) vs x and z (in units of a) inside the Z stack for 

I / Ic = 0.5 and b/a = 0.2. The contours correspond to magnetic field lines fiowing in 
the clockwise direction. The vertical dashed lines mark the boundaries of the middle 
region at \x\/a = c/a = 0.7292. 
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Figure 9. —Bz{x, z) in units of /ioa Jc/27r, calculated using a constant value of c at five 
heights in the stack, for a square stack {b/a = 1) with I /Ic — 0.5 and c/a = 0.6338a 
from (jB.5[) . The curves are offset vertically; by symmetry, Bz{0,z) = for any z, as 
marked by the horizontal dashed lines. 

Table 1. The constant c (in units of a) for I /Ic = 0.5 determined by the five different 
criteria discussed in the text. 

h/a i ii iii iv v 



oo 


0.5000 


0.5000 


0.5000 


0.5000 


0.5000 


10 


0.5286 


0.5163 


0.5163 


0.5218 


0.5081 


5 


0.5478 


0.5331 


0.5332 


0.5431 


0.5163 


2 


0.5895 


0.5803 


0.5816 


0.5927 


0.5416 


1 


0.6338 


0.6372 


0.6383 


0.6351 


0.5816 


0.5 


0.6809 


0.6996 


0.6821 


0.6727 


0.6383 


0.2 


0.7292 


0.7765 


0.7025 


0.6995 


0.6902 



the procedure given in Appendix B and ( IB. 51) . Let us call this criterion (i). We also 
could have determined c by choosing Bz{x, z) or B^xix, z) = dBz{x, z)/dx to be zero at 
various locations. For example, listed in Table [1] are values of c/a determined using the 
following criteria: (i) ([RS]), (ii) 5^(c,0) = 0, (iii) 5^:,.(0,0) = 0, (iv) 5^^(0,6/2) = 0, 
and (v) Bzx{S^,h) = 0. Table [2] gives values of the ratio 



for I/Ic = 0.5 and various values of b/a using the five different criteria. To show the 
sensitivity of Bz{x,0) to the choice of the constant c, we show in figure 10 plots of 
—Bz{x,0) and —Bz{x,b) vs x/a for b = a, I / Ic = 0.5, and values of c determined 
using three different criteria, (i), (ii) and (v). With criteria (i) and (ii), Bz{x,z) is 




(4) 
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Table 2. The average of Bz{x, z) over the quadrant 0<a;<ctoO<i;<6 divided by 
Bz{a, 0) for I/Ic = 0.5 determined by the five different criteria discussed in the text. 



b/a 


i 


ii 


iii 


iv 


V 


oo 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


10 


0.0000 


-0.0122 


-0.0122 


-0.0068 


-0.0204 


5 


0.0000 


-0.0149 


-0.0148 


-0.0047 


-0.0318 


2 


0.0000 


-0.0098 


-0.0084 


0.0034 


-0.0511 


1 


0.0000 


0.0039 


0.0051 


0.0015 


-0.0588 


0.5 


0.0000 


0.0215 


0.0014 


-0.0094 


-0.0486 


0.2 


0.0000 


0.0499 


-0.0275 


-0.0306 


-0.0401 



3 
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Figure 10. — ^^(a;, 0), —Bz{x,b/2), and —Bz{x,b) in units of for a square 

stack [h/a = 1) with I/Ic = 0.5, calculated using three different criteria (see text) to 
determine c/a: (i) 0.6338, (ii) 0.6372, and (v) 0.5816. For x/a < 0.5816, criterion ii 
yields the top curve, criterion i the middle curve, and criterion v the bottom curve. 
For x/a > 0.6372, the order is reversed, but all three criteria give very nearly the same 
value of Bz{x,z). The curves are offset vertically; by symmetry, Bz{0,z) = for any 
z, as marked by the horizontal dashed lines. 



very nearly equal to zero for |x| < c and z = but deviates more strongly sd, z = b, 
while the opposite is true for criterion (v). This provides additional evidence that the 
approximation of choosing c{z) to be independent of z can never yield Bz{x,z) — 
throughout the middle region, \x\ < c and \z\ < b. Note, however, that for c < |x| < a, 
Bz{x, z) is very nearly the same for all three criteria. 

For small b/a, the c = constant approximation is less successful in approximating 
the true fields in a Z stack. Figure 11 shows —Bz{x,0) along the center line {z = 0) 
for various currents with b/a — 0.5, and figure 12 shows the same for 6/ a = 0.2. In 
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Figure 11. —Bz{x,0) in units of /ioaJc/27r vs x/a for various currents in a stack of 
moderately low aspect ratio, b/a = 0.5. The values of c/a obtained from (jB.Sp are, for 
I/Ic = 0.2, 0.8889; 0.4, 0.7572; 0.6, 0.5952; 0.8, 0.3798; and 1, 0. (See figure 3.) 

the latter case the fields in the middle portion of the stack differ significantly from our 
desired condition = 0. Moreover, as can be seen in figure 12 for b/a = 0.2 and 
I / Ic = 0.2, calculations using criterion (i) yield values of —B^^x, 0) < for some values 
of X in the penetrated region c < x < a. Since these negative values are weighted by 
the factor (a — x) in ([2]), the losses calculated in the limit as ///c ^ using criterion (i) 
even become negative for b' = b/a < 0.0457, an unphysical result. On the other hand, 
we expect the errors in the losses due to such negative values of —Bz{x, 0) to be small 
for practical values of I/Ic > 0.2 when 6' = 6/a ~ 1. 

Figure 13 shows calculated values of —Bz{x, 0) vs x/a for various aspect ratios b/a. 
Note that for decreasing values of b/a, the profiles of ~Bz{x,0) become shallower and 
that the critical region where Jy = Jc is closer to the edges at x = a. It is clear from 
inspection of ([2]) that this behavior will lead to lower losses in Z stacks with smaller 
aspect ratios. We have used the above method for obtaining B^^x, z) in ([2]) to determine 
the hysteretic ac loss per cycle per unit length. 

Shown in figure [H] are our results for Q\ the loss per cycle per unit length, 
normalized to Q[^f (IHl), the loss per cycle per unit length of an infinite slab within 
a cross section Aab. Note that as the aspect ratio b/a increases, the ac loss per cycle 
converges slowly toward that of an infinite slab, but it can be significantly lower when 
6/a ~ 1. Included in figure [TH are our results for b/a = 0.2, although our approach 
becomes increasingly problematic at low aspect ratio, as discussed above. The solid and 
dashed curves show the dependence of the calculated loss on the criterion used to choose 
c/a. Clearly there is very little difference, which can be understood by inspecting the 
curves for —B^ vs x in figure 10. The behavior in the region |x| > c, which enters the 
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Figure 12. —Bz{x,0) in units of /ioaJc/27r vs x/a for various currents in a stack of 
low aspect ratio, b/a = 0.2. The values of c/a obtained from (jB.Sp are, for I/Ic = 0.2, 
0.9094; 0.4, 0.7967; 0.6, 0.6517; 0.8, 0.4489; and 1, 0. (See figure 3.) 




Figure 13. —Bz{x,0) in units of jioaJc/^Ti vs x/a for I/Ic = 0.5 and various aspect 
ratios, h/a = 0.2, 0.5, 1, 2, 5, 10, and oo (bottom to top). 
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Figure 14. Q' , the ac loss per cycle per unit length in a. Z stack in the anisotropic 
homogeneous- medium approximation, normalized to Q'l^^f the ac loss per cycle per 
unit length for an equivalent cross section of an infinite slab, using criterion (i) solid 
and (ii) dashed. The dotted curve shows the results of Norris 15J ([5]) for an isolated 
thin strip of thickness 2h = 0.4a. 



loss calculation, is only weakly influenced by the choice of c. 

For comparison, the dotted curve in figure \\M shows the Norris [15] result for Q'^j^jp, 
the hysteretic loss per cycle per unit length for an isolated thin, fiat strip of thickness 
26, normalized to Q'inf (ED) where 

Q'str^r> = ^^^^^^^[(1 - F) ln(l - F) + (1 + F) ln(l + F) - F'] (5) 
and F = I/I^. For F < 1, 

QWQW-(;^)f (6) 

The intercepts in figure 14 of the normalized ac loss in the limits I/Ic — > 1 and 
I / Ic — ^ are plotted vs b/a as the upper and lower curves in figure 15. The details 
of how to calculate Q' in these two limits are given in Appendixes C and D. When 
I/Ic Ij Q' is independent of the criterion used to determine the parameter c, and 
when I /Ic ^ 0, Q' depends only very weakly upon the criterion; in figure 15, the results 
calculated using criteria (i) and (ii) are almost indistinguishable. 

In the above calculations we have neglected the losses due to the currents and 
parallel ac fields {Br^) in the middle region |x| < c. In Appendix E we provide equations 
that can be used to estimate the middle-region losses. For typical tape dimensions and 
I/Ic > 0.2 we find that the middle- region losses are several orders of magnitude smaller 
than the losses at the edges (c < |x| < a). 




Figure 15. Q', the ac loss per cycle per unit length in a Z stack in the anisotropic 
homogeneous- medium approximation, normalized to Q'^^f ©, the ac loss per cycle per 
unit length for an equivalent cross section of an infinite slab, vs b/a. The upper curve 
shows the result Ri (|C.5P in the limit in the I/Ic 1, and the lower curves show the 
result i?o (|D-5|) in the hmit I/Ic ^ using criteria (i) soUd and (ii) dashed. 



3. Discussion and Summary 

Making use of what we have called the anisotropic homogeneous-medium approximation, 
we have introduced in this paper a theoretical framework for estimating the ac losses 
in a finite Z stack of superconducting tapes via straightforward analytic calculations. 
Our results yield Q', the hysteretic loss per cycle per unit length of a Z stack of total 
height 2b, where the tapes have width 2a. We have found it useful to compare our 
results with Q'^^f, the hysteretic ac loss per cycle per unit length for an equivalent cross 
section of an infinite slab. Our calculation is intended for application to the problem of 
calculating the hysteretic ac losses of a pancake coil, where the stack of tapes does not 
extend to infinity but rather curves back on itself. As long as the radius R of the coil 
is much greater than a or b, the field solutions should not differ significantly from what 
we have calculated, and the total hysteretic loss per cycle should be well approximated 
by 2tiRQ'. 

For a stack of many tapes, Q' is much larger than that for a single tape and 
this can be understood most simply by noting that Q' = '^Q'iniv where Q'init given 
in ([2]). Not only is the magnitude of the magnetic flux density Bz{x,z), which appears 
on the right-hand side of this equation, much larger than that generated by a single 
tape because of the superposition of the field contributions from all the tapes, but also 
Bz{x, z) is integrated in the z direction over a much greater height. For a stack of aspect 
ratio 6/a ~ 1 it is more appropriate to compare Q' with the hysteretic ac loss per 
cycle for an equivalent cross section of an infinite slab. For all finite values of b/a we 
find Q' < Q'inf- Referring to ([2]) and the field lines in Figs. 5-8, we see that the reason 
for this is that the field bends around the corners of the finite Z stack, and therefore 
\Bz{x, z)\ is always less than the corresponding quantity for the infinite slab. 

Grilli and Ashworth [2] have recently presented loss data for a multiturn pancake 
coil, which should be amenable to analysis using the formalism presented here. To 
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calculate the losses, they used a finite-element method, which apparently requires 
considerable computational resources. While we believe that our analytic approach 
has the advantage of efficiently yielding a loss prediction of sufficient accuracy for 
many applications, it would be of interest to compare results obtained using these two 
different approaches. One advantage of detailed finite-element calculations such as those 
discussed in |:2j is the capability of self-consistently incorporating the B dependence of 
the critical current density Jc{B). Our approach has made use of the assumption that Jc 
is a constant, independent of B. However, it would be possible for us to account crudely 
for the B dependence of Jc{B) by using a model for this dependence and replacing 
the constant Jc at each current amplitude / with Jc{Bm), where Bm is the maximum 
magnitude of Bz{x,z) at {x,z) = (a, 0). 

In section 2 and Appendix D we have pointed out that our anisotropic homogeneous- 
medium approximation is not accurate for small values of b/a, i.e., for a small number 
of superconducting layers in the stack. One reason for this is that we have assumed that 
each layer carries a constant average current density Jm in the middle region |x| < c. The 
motivation for this assumption is that B^ must be zero there and B^ between each pair 
of superconducting layers must be independent of x. However, the magnetic induction 
B^ at the top and bottom surfaces of the Z stack is not subject to this constraint but 
in general depends upon x, and consequently the screening sheet currents on the top 
surface of the top tape and bottom surface of the bottom tape also depend upon x. It 
seems likely that our failure to account for these sheet currents, which become relatively 
more important for a small number of layers, is the major reason for the problems with 
the present theory for small values of b/a. 
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Appendix A. A^, B^, and By for a cylinder of rectangular cross section 
carrying uniform current density 

In the following section we shall make use of several auxiliary functions. Consider the 
vector potential A{x,z) = yAy{x,z) and magnetic induction B{x,z) = xBx{x,z) + 
zBz{x,z) = V X A{x,z) generated by a uniform current density Jy in the region 
Xi < X < X2 and —b<z<b: 

Ay = {fioa^Jy/27i)fy{x[,x'2, b', x', z'), (A.l) 
B^ = {fioaJy/2TT)f^{x[, X2, b', x', z'), (A.2) 
B^ = {fioaJy/2']T)f;,{x[,X2, b', x', z'), (A. 3) 

where the primes denote dimensionless variables, x'l = Xi/a, etc. Since B^ = dAy/dx 
and Bz = —dAy/dz, we have fz = dfy/dx' and fx = —dfy/dz'. The dimensionless 
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functions fy, f^, and are (dropping the primes for simplicity): 

fy{xi,X2,b,x,z) = ^ - {x - xi){z + b)ln[{x - xif + {z + b^] 

+ {x - Xi){z - b) ln[(x - xif + {z- bf] 
+ (x - X2){z + b) ln[(a; - X2f + {z + bf] 

-{X - X2){Z - b) \n[{x - X2f + {z- bf] 

, .2 + \2 / ^ ~b \ 

— [x — Xi) arctan ( ] -\- (x — Xi) arctan f ) 

X — Xi^ X — Xi^ 

+ (x — X2) arctan ( ] — [x — X2) arctan ( ) 

X — X2' X — X2' 



— {z + bf arctan ( — ) + {z + b) arctan 



z + b ^ ' ' ^ z + b 

+{z - bf arctan ~ - {z - bf arctan (^— , (A.4) 

Z — b ' Z — n y-i 



fx{xi, X2, b, X, z) = {x — xi) In Y (x — Xif + {z + bf — {x — xi) In y (x — Xif + {z — bf 

— {x — X2) In \J{x — X2f + {z + bf + {x — X2) In \J {x — X2f + {z — bf 

+(2; + arctan ( — — + arctan ( — 

^ ' ^ z + b ^ ^ ' ^ z + b ^ 

— {z — b) arctan (- — + (z — b) arctan {- — , (A. 5) 

z — b ' z — b ' 

f,{xi, X2, b, X, z)^{z + b) In ^{x - xif + {z + bf -{z + b) In ^{x - X2f + {z + bf 



-{z-b)\n^{x- xif + {z- bf +{z-b)\n^{x- X2f + {z - bf 
+ [x — Xi) arctan ( ] — [x — Xi) arctan ( 

X — Xi'' X — Xi'' 

z + b \ z — b \ 
— ix — X2) arctan ( ] + [x — X2) arctan ( . (A. 6) 

X — X2^ X — X2'' 

Although these functions appear to have singularities in the In and arctan terms 
whenever the point {x,z) is on one of the boundaries {x = Xi, x = X2, z = —b, or 
z = b), the prefactors [{x — Xi), {x — X2), {z + b), or {z — b)] cause these terms to vanish 
there. Corresponding to the conditions that V • B = and J = V x B//xo, we have 
dfx/dx + dfz/dz — and 

dfx df. 



dz dx 



27r, Xi < X < X2 and \z\ < 6, (A. 7) 



= 0, otherwise. (A. 8) 

Appendix B. Finite Z stack for constant c 

We next wish to calculate the vector potential A(x, z) = yAy{x, z) and the corresponding 
magnetic induction 'Q{x,z) — xBx{x,z) + zBz{x,z) = V x A{x,z) generated by the 
following current densities Jy in a stack of height 2b: Jy — Jc ior c < \x\ < a and 
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Jy = Jm = 3m J c for |x| < c, wheie 

j„ = [l-(a/c)(l -///,)]. (B.l) 

The fields can be expressed as sums of contributions from tlie tliree regions, —a < x < 
— c, —c<x<c, and c < x < a, each of total height 2b, where —b<z<b: 

Ay{x, z) = (/ioaVc/27r)[/y(-l, -c', 6', x', /) + jmfy{-c', c', 6', x', z') + fy{c', 1, b\ x\ z!)\ 

(B.2) 

B^{x,z) = (/xoa-^c/27r)[/^(-l, -c' ,b' , x' , z') + jmfx{-c' , c' ,b' , x' , z') + fx{c', l,b' , x' , z')], 

(B.3) 

B,{x, z) = (/ioaJ,/27r)[/,(-l, -c', 6', x', + j„/.(-c', c', 6', x', z') + /,(c', 1, 6', x' , z% 

(B.4) 

where x' = x/a, = 2;/a, b' = b/a and c' = c/a. For given values of b' = b/a and ///c, 
we determine c using criterion (i) by requiring that 

pb re rb 

dz dxB^{x,z)= / dz[Ay{c,z)-AJO,z)]=0, (B.5) 
JO Jo Jo ' 

which can be solved numerically, using analytic expressions for the integrals of fy 
required in ( IB. 51) . Figure [3] shows plots of c' = c/a, obtained from ( IB.ip . ( IB. 21) . and 
(]B.5|) . vs I/Ic for various values of b' = b/a, and figure H] shows corresponding plots of 

Jm- 

Contours of constant Ay{x,z) obtained from flB.2p correspond to magnetic field 
lines. 



Appendix C. Losses for I/Ic 1 

In the limit as z = ///c ^ 1 [see figure 3], c independent of either the value of 
b/a or the criterion (i)-(v) used to determine c. We can then apply the anisotropic 
homogeneous-medium approximation to calculate Q' = '^Q'init from (2) and flB.4p using 
c' = and /^(O, 0, b' , x' , z') = 0. The result for I/I^ = 1 is 



where 



Q' = -fioJXfib/a), (C.l) 

71 

f{u) = ^ [ - 3m^ + 8m(1 - u^) tan"^(M) - 2m(3 - Au^) tan"^(2M) 



Expansions of f{u) about u = oo and u = yield the leading terms, 

f{u) = nu/3, u ^ oo, (C.3) 

= 2(ln4-l)M^ u^O. (C.4) 
When I/Ic = 1, we obtain from (IC.ll) and (3) the ratio 

. «^ . ,c,, 
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which is plotted vs b/a as the upper curve in figure 15. When I/If. = 1 and b/a ^ oo, 
we find from flC.3l) and flC.51) that Q'/Q[^j 1, as expected. When I/Ic = 1 and 
b/a 0, we find from (EH), (El), and (ED that Q'/Q'^trip ^ 1- 

Appendix D. Losses for I / Ic ^ 

According to the above anisotropic homogeneous-medium approximation with c 
independent of z, the ac loss per cycle per unit length Q' = ^Q'init (E]) is proportional 
to {I / Ic)^ in the limit as i = I /Ic — ^ 0. In this limit we have c' = 1 — e and [from ( IB. II) 
Jm = — e to first order in e, where e = ki and k is a constant of order unity {k < 1), 
which depends upon the criterion used to determine c. By expanding fy in (IB. 21) or fz 
in (IB.4P to first order in e and carrying out the integration in (IB.Sp . we obtain for the 
value of k using criterion (i) 

ki = [87r6'3-486'2ln2- 166'(3-6'2)tan-H' + 86'(3-46'2)tan-i(26') 
+ 8(l-36'2)ln(l + 6'2) - 2(l-126'2)ln(l+46'2)]/ 
[87r6'3 + 166'3tan-i b'- 8b'{3+Ab''')tan-\2b') + 126'2ln(6'2) 
- 4(l + 36'2)ln(l + + 41n(l+46'2)]. (D.l) 

Similarly, by expanding fz in (IB. 41) to first order in e and setting -Bz(c, 0) = 0, we obtain 
for criterion (ii) 

4tan-i(&72) + yin(l+4/y^) 
~ 7r + 2tan-i(&'/2) + &'ln(l+4/6'2)" ^ ' 

Carrying out the integration required in ([2]), we obtain in the limit as i = ///c — > 0, 

Q' = -fioJyg{k,b'){I/hf, (D.3) 

71 

where 

g{k,b') = [46'tan-H' + 6'2ln(l + l/6'')-ln(l + 6'2)]A;2 

- [nb'/3 + 26' tan-i 6' + 6'" ln(l + 1 /fo'^)] fcl (D.4) 
The ratio of the result in (ID.3P to Q-^j ([3]) is 



^ - " 

which is plotted as the lower curves in figure 15 for two of the criteria [(i) solid and (ii) 
dashed] used to determine the constant c. As b' ^ oo, ki ^ 1, ka ^ 1, and Rq ^ 1 for 
both criteria, and as b' ^ 0, k^ ^ 0, ka 0, and Rq for both criteria. However, 
for very small values of b', where the present approach is not accurate, the value of Rq 
for k = ki is negative for < 6' < 0.0457, an unphysical result. 
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Appendix E. Contribution to the losses from the middle region, |x| < c 

In the above sections we have calculated the ac losses in the outer regions (c < |a;| < a) 
due to perpendicular magnetic flux (oc B^) moving in and out from the edges of the 
tapes. We have so far neglected the ac losses in the middle region < c) due to 
parallel magnetic flux (oc B^) moving in and out from the top and bottom surfaces 
of the tapes, on the assumption that these losses are very small. In this appendix we 
present equations that can be used to conflrm that the losses in the middle region are 
indeed much smaller than those in the outer regions. 

Consider the flnite Z stack sketched in flgure 1, and label the z coordinate of a 
given tape as z„. For an odd number of tapes in the stack Zn = nD, n = 0, ±1, ±2, 
and for an even number of tapes Zn = ±D/2, ±3D/2, ±5D/2,... . In the middle 
region < c, each tape carries an ac current of reduced amplitude jm = Jm/Jc, 
normalized to the critical current density. However, since dB^/dz = fioJm, we have 
Bx = fioJmZ, such that the tape at z„ is also subject to an applied in-phase ac magnetic 
induction of amplitude /ioJmknl- It is convenient to normalize this to the penetration 
fleld Bp = fiojcd/2 = fiQJcD/2, such that the applied fleld has the reduced amplitude 

hn = im\Zn\/{D/2). 

Q^, the hysteretic loss per cycle per unit volume of a superconducting slab subjected 
to an ac parallel fleld of reduced amplitude h and an in-phase current of reduced 
amplitude i = j/jc < 1, has been calculated by Carr [18], whose results can be expressed 
as 

Qv = ^f{h,z), (E.l) 



where 



3/xo 



f{h,i) = i^{i + 3h), h < i, (E.2) 
= h{h'^ + 3i'^), i<h<l, (E.3) 

= /i(3 + - 2(1 - z^) + — ^ -J-,h>l. (E.4) 



(h-i) {h-if 

q'^, the hysteretic loss in tape n per cycle per unit length due to the ac current and 
parallel fleld in the middle region (|x| < c) of the stack, can therefore be calculated from 

AB'^cd 

q'^ = Q,2cd = —I^f{h^,j^), (E.5) 

3yU0 

and Qm, the total hysteretic loss per cycle per unit length of the Z stack due to the ac 
currents and parallel flelds in the middle region (|a;| < c), can be obtained by carrying 
out the sum over all tapes, 

Q'm = T.C (E.6) 

n 

As can be seen from (IE.2p - (lE.4p . the appropriate expression for f{hn,jm) to be used in 
the sum depends upon the value of hn relative to jm and 1. The ratio of Q'^ to Q'^^f & 
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Figure El. R,no (E.8), the contribution Q'^^ (E.6) to the ac loss per cycle per unit 
length due to the middle region {\x\ < c) of the Z stack, normalized to ([3]), the 

ac loss per cycle per unit length for an equivalent cross section of an infinite slab, in 
the limit as I/Ic 0, vs N, the number of tapes in the stack, for d' = 2 x 10~^, 
D' = 2 X 10^^, and N ~ 1006' ~ 100b/ a. The results were calculated using two 
different criteria [(i) solid and (ii) dashed] for calculating the parameter c. 



IS 



Numerical evaluation shows that, in contrast to the behavior of Rq and -Ri shown in 
figure 15, Rm is a monotonically decreasing function of I/Ic with its maximum value at 
I/Ic = 0, given (for an even number N of tapes) by 

R^o = d'D'N{N^ + 4) (1 - fc) Vl6, (E.8) 

where d' = d/ a, D' = D / a, and k [see (ID. II) and (ID. 21) ] depends upon the criterion used 
to determine c. To estimate the order of magnitude of the middle-region losses we use 
the following assumptions: 2a = 10 mm, D = 100 /im, d = 1 fim, 2b = ND, such that 
d' = 2x 10-^ D' = 2x 10"^ and = 1006' = 1006/a. Shown in figure [El] is a plot 
of Rmo VS N, the number of tapes in the stack, for the two criteria [(i) solid and (ii) 
dashed] we have used to calculate c. Figure [E2l calculated for b' = 1, shows the general 
behavior of how the ratio Rm/Rmo depends upon I/Ic- This ratio has its maximum as 
I/Ic 0, where c' = c/a ^ 1, but vanishes as I/Ic 1, where c' = c/a — 0. Although 
both criteria (i) and (ii) for choosing c were used to calculate this ratio, the two curves 
are indistinguishable on this plot. Comparing figures lETl and [E2] for the middle-region 
losses with figures 14 and 15 for the outer-region losses, we see that for I/Ic > 0.2 
the hysteretic losses from the middle region of the stack |a;| < c are typically at least 
two orders of magnitude smaller than those from the edges of the tapes. Howeover, in 
the limit as I/Ic 0, where the middle region includes nearly the entire volume and 
the outer regions shrink to zero, the middle-region losses become more important but 
still remain relatively small. In summary, these results confirm that for I/Ic > 0.2 the 
hysteretic losses from the middle region of the stack |x| < c are typically several orders 
of magnitude smaller than those from the edges of the tapes. 
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Figure E2. The ratio Rm/RmO, calculated from (E.7) and (E.8), vs I/Ic assuming 
rf' = 2 X 10"^, D' = 2x 10-^ and N = 100 or b' = b/a = 1. Two different criteria [(i) 
and (ii)] for calculating the parameter c yield nearly identical results. 
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